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Abstract 



This work concerns superharmonic perturbations of a Gaussian measure given by a special class of positive 
weights in the complex plane of the form w(z) — exp(— \z\ 2 + U^(z)), where U^{z) is the logarithmic potential 
of a compactly supported positive measure p. The equilibrium measure of the corresponding weighted energy 
problem is shown to be supported on subharmonic generalized quadrature domains for a large class of perturbing 
potentials t/ M (z). It is also shown that the 2x2 matrix d-bar problem for orthogonal polynomials with respect to 
such weights is well-defined and has a unique solution given explicitly by Cauchy transforms. Numerical evidence 
is presented supporting a conjectured relation between the asymptotic distribution of the zeroes of the orthogonal 
polynomials in a semi-classical scaling limit and the Schwarz function of the curve bounding the support of the 
equilibrium measure, extending the previously studied case of harmonic polynomial perturbations with weights 
w(z) supported on a compact domain. 



1 Introduction 

This work mainly concerns equilibrium problems in potential theory, but its motivation derives largely 
from two related domains: random matrix theory and interface dynamics of incompressible fluids. Recent 
work of P. Wiegmann, A. Zabrodin and their collaborators [18, 20] connected the spectral distributions 
of random normal matrices to the Laplacian Growth model for the interface dynamics of a pair of two- 
dimensional incompressible fluids. The unitarily invariant probability measure on the set ofnxn complex 
normal matrices in [18] is determined by a potential function 

V(z,z) = zz + A(z) + Ajz) (1.1) 
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and the corresponding density is of the form 

f(M > M*) = e-* TrV ( M > M *) (1.2) 

with respect to a unitarily invariant reference measure on the set of normal matrices. The function A(z) 
is assumed to have a single-valued derivative A'(z) meromorphic in some domain D C C. Under suitable 
assumptions on A(z), the large n limit of the averaged normalized eigenvalue distribution in the scaling 
limit 

n — > oo, H — ► 0, nh = t (1-3) 

tends to a probability measure fiy,t where t is some fixed positive number {quantum area). It turns out 
that fiyt is the unique solution of a two-dimensional electrostatic equilibrium problem in the presence of 
the external potential V in the complex plane. In most cases, nv,t is absolutely continuous with respect 
to the planar Lebesgue measure with constant density. It can be shown that the support supp(/iyt) 
undergoes Laplacian growth in terms of the scaling constant t and the area of the support is linear in t. 
However, this evolution problem is ill-defined; even initial domains with analytic boundaries may develop 
cusp-like singularities in finite time and the solution cannot be continued in the strong sense (see [18], 
[17]). The averaged eigenvalue density of the corresponding normal matrix model for finite matrix size 
n be viewed as describing a sort of discretized version whose contiuuum limit may be interpreted as a 



semiclassical limit (1.3 1 that tends to /iy,t as shown in [1,9]. 

In studying these questions, it is important to understand first the possible shapes of compact sets 
which are supports of equilibrium measures for potentials of the form ( |1.1[ ). In this work we show the 
support of the equilibrium measure for a class of perturbed Gaussian potentials of the form 

V a , u (z):=a\z\ 2 + U"{z), (1.4) 

where a > and v is a compactly supported finite positive Borel measure, to be so-called generalized 
quadrature domains [8,14]. 

To understand the asymptotic behaviour of the averaged eigenvalue density of normal matrix models 
in different scaling regimes, one has to study the asymptotics of the corresponding orthogonal polynomials 
for the weight e~ v ^ z \ The well-known Riemann-Hilbert approach is not applicable directly in this case 
because the orthogonality weight is not constrained to the real axis. However, there is a sort of matrix 
<9-bar problem, introduced by Its and Takhtajan [10], which is a candidate to replace the matrix Riemann- 
Hilbert problem in the study of strong asymptotics of the orthogonal polynomials. In the present work 
in this regard is shown that the 9-bar problem is also well-defined and characterizes the orthogonal 
polynomials for the class of perturbed Gaussian potentials considered above. 

To fix notations, let m denote the two-dimensional Lebesgue measure in the complex plane C. We 
denote respectively by H and H c the closure and the complement of a set H C C and by Ig the indicator 
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function.. The open disk of radius r centered at c € C is denoted by B(c, r), and the Ricmann sphere by 
C. 



2 Weighted energy problem and logarithmic Potentials 

In this section we briefly describe both the classical and weighted energy problems of logarithmic potential 
theory (see [12] and [15]) and specify the class of background potentials we are concerned with in this 
paper. 

Definition 2.1. Let \i be a compactly supported finite positive Borel measure in the complex plane. The 
logarithmic potential produced by /x is defined as 

U»{z) := /"log— !— rd/iM (zeC). (2.1) 
Jc \ z ~ w \ 

In particular, for a bounded subset S of the plane with m(S) > 0, we consider the measure T]s given 

by 

dvs = I s dm. (2.2) 

Thus 77s is the Lebesgue measure restricted to S. In the following, we use the simplified notation II s (z) 
for the logarithmic potential U vs (z) of the measure r/s- 

The logarithmic potential of a positive measure \x is harmonic outside the support of /z and superhar- 
monic on supp(/j) (see [15], Theorem 0.5.6). Moreover, it has the asymptotic behaviour 

U»(z) log^ + o(-), (M-oo) (2.3) 



\z\ \z 

where /x(C) is the total mass of \i. If U^{z) is smooth enough the density of the measure /a can be 
recovered from this potential by taking the Laplacian of Z7 M (z): 

Theorem 2.2 ( [15], II. 1.3). If in a region RCC the logarithmic potential U fl (z) of the measure \x has 
continuous second partial derivatives, then [i is absolutely continuous with respect to the planar Lebesgue 
measure m in R and we have the formula 

dn= -^-AU^dm. (2.4) 

Now, let K be a compact subset of C and let M. (K) denote the set of all Borel probability measures 
supported on K. In classical potential theory, the logarithmic energy of a measure \i € A4(K) is defined 
to be 

/( M ):= / U»{z)dn{z)= ( ( log— ^—d^d^z). (2.5) 

The quantity 

Ek:= inf (2.6) 
neM(K) 
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is either finite or +00. The logarithmic capacity of K is 

cap(K) := e~ Elc . (2.7) 

If Ek < 00 then, by a well-known theorem of Frostman (see e.g. [12]), there exists a unique measure \±k 
in M.(K) minimizing the energy functional /(•) and this measure is called the equilibrium measure of K. 
The capacity of an arbitrary Borel set B C C is defined as 

cap(B) := sup{cap(i^) | K C B,K compact}. 

A property is said to hold quasi- everywhere if the set of exceptional points (i.e. those where it does not 
hold) is of capacity zero. 

In the more general setting we have a closed set ECC and a function w : £ — > [0, 00) on £ called the 
weight function. Usually the weight function is given in the form 

w(z) = exp(-Q(z)) (2.8) 

where Q: £ — > (—00,00]. In the electrostatical interpretation the set £ is called the conductor and Q is 
called the background potential. 

Definition 2.3 ( [15]). The weight function w is said to be admissible if 

• w is upper semi- continuous, 

• {z G £ I w(z) > 0} has nonzero capacity, 

• lim |z|u>(z) = 0. 

\z\ — >oo 

The admissibility conditions can be rephrased in terms of the potential Q; w(z) — exp(— Q(z)) is 
admissible if and only if Q is lower semi-continuous, the set {z G £ | Q(z) < 00} has nonzero capacity 
and lim (Q(z) — log \z\) = 00. 

\z\— >oo 

Let Al(£) denote the set of all Borel probability measures supported on £ C C. The weighted energy 
functional Iq is defined for all \i e M(T<) by 

I Q (H) := JJlogftz-tWzMt^dnWdnit) (2.9) 

= JJjogj^-jd f i(z)d f i(t) + 2 J^Q{z)d l i{z). (2.10) 

The goal is then to find a probability measure that minimizes this functional on Al(£). If Q is admissible 
it can be shown (see [15], Theorem 1.1.3) that 

E Q := inf ; /q(m) (2.11) 
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is finite and there exists a unique measure, denoted by /uq, that has finite logarithmic energy and min- 
imizes Iq. Moreover, the support of /iq, denoted by Sq, is compact and has positive capacity. The 
measure \iq is called the equilibrium measure of the background potential Q. The logarithmic potential 
satisfies the equilibrium conditions 

U tlQ {z) + Q(z) > Fq quasi-everywhere on S, 

U^{z) + Q{z) < F Q for all z e Sq, (2.12) 
where Fq is the modified Robin constant: 

Fq = E q - j QdiiQ. (2.13) 

Motivated by random normal matrix models (see [18], [17]), we are interested in background potentials 
of the following type: 

V a , v (z) :=a\z\ 2 + U v {z), (2.14) 

where a is a positive real number and v is a compactly supported finite positive Borel measure. These po- 
tentials have a planar Gaussian leading term controlled by the positive parameter a and this is perturbed 
by a fixed positive charge distribution given by the measure v. 

Proposition 2.4. The potential V a , u (z) is admissible for all possible choices of a and v. 

Proof. V a . u (z) is lower semi-continuous because U v (z) is superharmonic in the whole complex plane. 
The set where V a ^(z) is finite contains at least C \ supp(^), which is of positive capacity since supp(^) 
is compact. Finally, the required boundary condition is also fulfilled: 

V a , v (z)-iog\z\ = a|z| 2 + [r(z)-log|z| (2.15) 
= a\z\ 2 -(v(C) + I) log \z\ + o(^j , (2.16) 

so the difference V a ^(z) — log \z\ goes to +00 as \z\ — > 00. 

□ 

We are especially interested in cases for which the perturbing measure v is singular with respect to 
the planar Lebesgue measure to. In particular, v can be chosen to be a positive linear combination of 
point masses, i.e. 

m 

^:=]T/3 fe <Sa fc , (2.17) 
fe=i 

where a l7 a 2 , . . . , a m £ C are the locations and 0i, fei ■ ■ ■ > Pm are the charges of the fixed point masses 
respectively. 



5 



3 Supports of equilibrium measures and quadrature domains 

The determination of the support of the equilibrium measure for a background potential V(z) — V a _ v {z) 



of the form (2.141 above is closely related to finding generalized quadrature domains of some measures in 



the complex plane. Let us recall the definition of quadrature domains given by M. Sakai ( [14]). 

Definition 3.1. Let v be a positive Borel measure on the complex plane. For a nonempty domain (open, 
connected) CI in C let F(Cl) be a subset of the space 

ReL 1 ^) := {Re/ | / £ L 1 (Q, 1 m)) . (3.1) 

of real-valued integrable functions on CI. 

The domain CI is called a (generalized) quadrature domain of the measure v for the function class 
F(Cl) if 

(a) v is concentrated on CI, i.e. v{Cl c ) = 0. 
(b) 

/ f+dv < oo and [ fdv < [ fdm (3.2) 
Jn Jn Jn 

for every f £ F(Cl) where f + :— max{/, 0}. 

Note that if F(Ct) is a function class such that — / £ F(Ct) whenever / £ F(£l) then the second 
condition is equivalent to 

/ \f\dv < oo and / fdv = I fdm (3-3) 

Jn Jn Jn 

for every / £ F(f2). We are interested in the following subclasses: 

ReAL^fl) = {Re/ £ L 1 ^,™) \ f is holomorphic in Q} (3.4) 
HL 1 ^) = {h £ L 1 ^,™) | h is harmonic in fl} (3.5) 
SL 1 ^) = {s £ L 1 (fi,m) | s is subharmonic in ft} (3-6) 

For a measure v the quadrature domains corresponding to these classes are called classical (holomorphic), 
harmonic and subharmonic quadrature domains respectively. We have the obvious inclusions 

Q(u, SL 1 ) C Q(v, HL 1 ) C Q(v, ReAL 1 ), (3.7) 

where Q{v, F) denotes the set of quadrature domains of v for the function class F . It is important to 
note that if the domain f2 belongs to Q(y, F) then its saturated set or areal maximal set 

[CI] :={z£C\m (B(z, r) n Cl c ) = for some r > 0} (3.8) 



G 



also belongs to Q(i/,F). 

For example, it can be shown that the disk B(c,R) is the only classical quadrature domain for the 
point measure v = R 2 ir5 c (see [14], Example 1.1). The simplest examples are the classical quadrature 
domains whose quadrature measure is a positive linear combination of point masses: 

m 

f = P k eR + . (3.9) 

k=l 

This means that for every holomorphic function / that is integrable on we have the identity 



« m 
/ fdm = y^(3 k f(a k ). 
Jn k=i 



(3.10) 



An immediate generalization is obtained by allowing points a k of higher multiplicities fk > 1 in the above 
sum, this means allowing derivatives of finite order to appear in the sum representing the area integral 
functional: 

/. m jk 

/ fdm = J2J2^jf W (a k ). (3.11) 
^ n k=l 1=0 

However, in this work we do not consider such quadrature domains. 

It is easy to see that if Q is a subharmonic quadrature domain of the measure v then, using subharmonic 
test functions of the form 

«*(«;)= -log— !— r (zGC), (3.12) 
\z — w\ 



we have 



U n {z) < U v (z)iizeC, 

U n (z) = U"(z) if zeC\n. (3.13) 



To illustrate the structure of the equilibrium measure of a potential of the form (2.141 above, we consider 
a simple but nontrivial example: 

V(z) = a|z| 2 + /31og^^, (3.14) 

where a £ , £ R + and a £ C. The calculation of the equilibrium measure for this potential is quite 
standard (see, for example, [18], [17]) but the details will be of importance in suggesting generalizations. 
For the sake of completeness, the statement of the result and a short sketch of its proof are therefore 
included here. 

To find the equilibrium measure for this potential one can use the following characterization theorem: 

Theorem 3.2 ( [15], 1.3.3). Let Q: £ — * (—00,00] be an admissible background potential. If a measure 
a £ has compact support and finite logarithmic energy, and there is a constant F £ K such that 

U a (z) + Q(z) = F quasi- everywhere on supp(o~) (3.15) 
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and 

U a (z) + Q(z) > F quasi- everywhere on E, (3.16) 
then a coincides with the equilibrium measure \iq. 

The logarithmic potential of the uniform measure rj B ^ c R j on a disk B(c, R) is easily calculated to be 



U B{c > K) {z) = I (3.17) 



.R 2 7rlogi r |z-c|>it 

Z - C 



Proposition 3.3. Define two radii R and r as 



R:—\ and r := \ — . (3.18) 

V 2a V 2a v ' 

The equilibrium measure /iy is absolutely continuous with respect to the Lebesgue measure and its density 
is the constant . The support Sy of \xy depends on the geometric arrangement of the disks B(a,r) and 
B(0, R) in the following way: 

(a) If B(a,r) C 5(0, R) then 



S v = B(Q,R)\B(a,r) (3.19) 



(see (a.l) and (a.2) in Figure 3.1l. 
(b) If B(a, r) B(0, R) then C \ Sv is given by a rational exterior conformal mapping of the form 

/:C\{C : |C|<l}^C\Sy, /( C ) = p C + u+-^, (3.20) 
where the coefficients p € K + ,0 < \A\ < 1 andu,v G C of the mapping /(C) are uniquely determined 



by the parameters a, (3 and a of the potential V(z) (see (b.l) and (b.2) in Figure 3.1l. 



Proof. Suppose first that B(a,r) C B(0,R). 
Let a be the measure given by 



da := -^—I K dm (3.21) 
m(K) 



where K = B(0, R) \ B(a, r). The area of K is 



m(K) = (R 2 -r 2 )n= — . (3.22) 
2a 



Therefore the logarithmic potential of a is 

2a 



(z) = — (u B ^ R \z) - U B <* r \zj) . (3.23) 



8 





(b.l) 



(b.2) 



Figure 3.1: The shape of the support Sy (shaded area) is illustrated for the different disk configurations 
(a) and (b) in Proposition 3.3 



Now, for z E K the effective potential at z is 



U°{z) + V{z) = «i? 2 log — + 1-U_ -2ar 2 log 



R 2 R 2 



\z-a\ 



+a\z\ 2 + /?log 



\z — a\ 



1 



= aRf ( log + 1 



(3.24) 



Define 



F~nR- (log-^ + 1 



(3.25) 
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If z $ K a short calculation gives 

U a (z) + V(z) = 

where 



F + 2ar 2 f 
F + 2aR 2 f 



z e B(a, r) 
\z\ > R. 



x 2 -l 



- logs. 



(3.26) 



(3.27) 



Since / is nonnegative the effective potential satisfies 

U a (z) + V(z) > F. (3.28) 

By Theorem |3.2[ we conclude that the equilibrium measure for the background potential V is a. 

Now suppose that B(a,r) (£ B(0,R). For this case, we only sketch the calculation giving the system 
of equations that relates the parameters of the potential V(z) and the conformal map /(C)- The potential 



is smooth in the domain C\ {a} and its Laplacian there is constant which suggests, by Theorem 2.2 that 
the density of the equilibrium measure is equal to 



(3.29) 



We expect the equilibrium measure therefore to be the normalized Lebesgue measure restricted to some 
compact set S c C \ {a}; that is, 



1 



m(S) 



dm 



where m(S) is given by (3.22). The first equilibrium condition (3.151 is 

U llv (z) + V(z) = F (zeS) 



(3.30) 



(3.31) 



for some constant F. Assuming the necessary smoothness of (z) and applying the differential operator 
d z , this gives the necessary condition 

a f dm(w) 

7T 



_ (3 1 

az 



which is the same as 



a f dw A dw 



2iri 



By Green's Theorem we get 



1 

2wi 



wdw 
,iS 7 -- w 



2 z — a 
2 z — a 

2a z — a 







(z G S), 
= (zeS). 

(z G S). 



(3.32) 
(3.33) 

(3.34) 



Following [18], we seek to express the compact set S in terms of an exterior conformal mapping of the 
form 



/:C\{C : |C|<1}->C\5, f(()=p( + u + 



C-A 



(3.35) 
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where p > and < \A\ < 1. By rewriting the previous form of the derived equilibrium condition we 
obtain 

i f /(C)f(CR_ P i {zesy (3 . 36) 



2ni J\(\=i z — /(C) 2a z — a 

Along the positively oriented simple circular contour |£| = 1 we have 



/(C)=P^ + «+ I ^. (3.37) 



Therefore the rational function 

z-pC-u- - C _ A 
has coefficients depending on p, m, u, A and satisfies the equation 



T(C; z) := ^ 7" s/ v „ (C A)2/ (3.38) 



^ [ T(C,z)d( = -f— (zeS) (3.39) 
2?r* i|c|=i zaz-a 



subject to the area normalization condition (3.221. 

The differential T((; z)d( has four fixed poles in £ at £ = 0, £ = oo, ( = A, ( = = and two other 
poles depending on z via the equation z — /(C)- (These extra poles may coincide with some of the fixed 
poles above.) This equation can be rewritten as a quadratic equation in 

p( 2 + (u~z- Ap)( + A(z-u)+v = 0. (3.40) 

If z £ S both solutions of this equation are inside the unit disk {£ | |£| < 1}. To calculate the integral of 
T(£; z)d( in terms of residues, we write the the contour integral in the standard local coordinate £ = | 
around £ = oo: 

If If ^{1 \ d£ 



., ■ , T{Q;z)dQ=— / T -|. (3.41) 

2?" J|C| = 1 27TZ 7| ? | =1 V? / ? 

The two simple poles inside the disk {£ | |£| < 1} are £ = and £ = ^4 and hence 

^jL t (H P =Res ^° T (H ^ +Res ^ r (^) f ^ 

Since 



the residues are 



r( i U + (p-(Af)(i-^) 



nes._ n r(i;*j 1 = =^-7/. (>!.-]-]) 



1 \ 1 



(/'-fwf)(iHA| 2 ) 



Res, -tT — = = ^ — L - L ^ ~. (3.45) 

*- A [ Cje A { i^\A\>){zA-uA-p)-vA 2 ^ ' 
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A short calculation gives the area of S in terms of the mapping parameters: 

m {S)=*( P 2 - n lV }' )- (3-46) 



(1-L4| 2 ) 2 

Finally we obtain the following system of equations: 



■2 



1 



. i /1 1 2 

-u = 



(_1 - |A| 2 ) 2 2a 
v 



A 

p vA (3.47) 

u + = + -— = a 

A l-\A\ 2 

v ( vA 2 \ (3 

We must prove that if we assume \a\ + r > R there exists a unique solution p, u,v, A to this system in 
terms of the parameters a, (3, a. Eliminating u from the third equation gives 

" = 3 + A(T^W < 3 ' 48 > 



The last equation of (3.471 shows that -4y is real. The phases of v and A are therefore fixed by (3.48). 
Writing a in polar form 

a = te^ (3.49) 
(with t > because B(a, r) <£_ B(0, R)), we obtain 

w^se 2 ^, A = Ke i<t ', (3.50) 
where s and K are positive real numbers. We can express p and s in terms of K: 

P = ^(^ + 1) (3.51) 



2Kt V 2a, 



Setting x = K , this must be a solution of the the cubic equation 



2tV - ( t 4 + ^^i 2 ") x 2 + i - 0. (3.53) 



The condition lal + r > i? means that 



+4 1 + 2/3 2 , 1 



r -t 1 + — < 0. (3.54) 

a 4a 2 



Defining the function 



:=2/- l ,:' i - + J_ (3.55) 
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wc have 



5 (0) = > and g(l) = t A - 1 ± 2 l t 2 + < (3.56) 
4a- 2 a 4a z 



</(x) = 6t 4 x 2 -2[t 4 + ] -^-t 2 

a 



.r 



by (3.541. Since 



= ^-^r+^JJ (3 ' 57) 

is negative in the interval [0, 1], g(x) has a unique root in (0, 1), and herefore K is uniquely determined 



by (3.531. This means that there is a unique solution for p,u,v and A of (3.47) in terms of a,f3, a. 

To conclude the proof one should show that the logarithmic potential — II s (z) satisfies the inequality 



(3.161 of Theorem 3.2 This part of the proof is omitted. 

□ 



The electrostatic interpretation of case (a) in Proposition 3.3 is simple. If we replace the point charge 
(35 a by a uniform charge distribution of density — on B(a, r), the resulting configuration is in equilibrium 
in the presence of the pure Gaussian potential a\z\ 2 . The disk B(a,r) is a quadrature domain for the 
measure r 2 TrS a , so 

— U B(a ' r) (z) = /3U S °(z) in Z £%r) c , (3.58) 

7T 

which means that the electric fields of ^?7_B( a .r) and (3S a are indistinguishable in the exterior of B(a, r). A 
quadrature domain shaped cavity emerges in the support of the equilibrium measure of the unperturbed 
Gaussian potential since the fixed perturbing measure substitutes the portion of uniform charge placed 



in the cavity of the original equilibrium configuration, as illustrated in (a.l) and (a.2) of Figure 3.1 A 
useful generalization of this idea turns out to be valid in a more general setting: 

Theorem 3.4. Let 

V(z) :=a\z\ 2 +U v (z) (3.59) 
be a background potential. Assume that the measure v can be decomposed into a sum 



m 

V = 

k=l 



5>fc (3.60) 



where the measures v\. are all finite positive Borel measures satisfying the following conditions: 

(a) The supports of the measures Vk are pairwise disjoint and each v% has positive total mass. 

(b) The measure -^Vk has an essentially unique (i.e. unique up to sets of measure zero) subharmonic 
quadrature domain, and D). denotes the saturated element of Q(uk, SL 1 ) for all k = 1,2, ... ,m 
respectively. 
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(c) The domains Dk are pairwise disjoint and D k C B(0, R) for all k = 1,2, ... ,m where 

R := 



1 + W 



2a 



(3.61) 



Then the equilibrium measure [ly is absolutely continuous with respect to the Lebesgue measure with 
constant density — and is supported on the set 



K:=B(0,R)\ (J D k 



(3.62) 



\k = l 



(The situation is illustrated for a simple configuration of point and line charges in Figure 3.2 ) 




Figure 3.2: A typical configuration involving subharmonic quadrature domains: a disk, a so-called bicir- 
cular quartic (see [16]) and an ellipse 



Proof. Let a be the measure given by 



da := 



1 



-Ixdm 



m(K) 

To calculate the area of K, we note that the area of Dk is given by 



m{D k ) 



dm = — / dvu — — za,(C). 
Dk 2a J 2a ' 



Therefore 



m(K) = m(B(Q,R))-^m(D k ) 



k=l 



Mi + KQ-X>(c) 



2a 
2a' 



k=l 



(3.63) 
(3.64) 

(3.65) 

(3.66) 
(3.67) 
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For each measure Vk the corresponding logarithmic potential XJ Vk (z) satisfies 

2a 



U Dk (z) < U" k (z) if z E C, (3 
— U Dk (z) = U Vk (z) if z£C\D k . (3 

7T 



The logarithmic potential of a is 



U°(z)= — (u B ^ R \z)-f^U D *{z)\ . (3 



\ fe=i 
Now, if z e if then the effective potential at z is 



cr(*) + v[z) = ^(log^ + i-^-^E^M-) 

+a|z| 2 + £T(z) 



aR 2 (log i + 1 ) - £ + U»(z) 



fe=i 



= «i? 2 (log ~ + 1 J . (3 

Define F := aR 2 (l og -gj + l). If z (j[ K then either z E for some k or |z| > R. If z G then 
have the inequality 

U°{z) + V{z) = a R*L g ^ + l-^)-^jru Dk (z) 

+a\z\ 2 + U v {z) 
= F+C/' yfc (z)- — U Dk (z) 

7T 

> F. (3 ' 

On the other hand, if j jar j > R then 

1 9 rn 

U a (z) + V(z) = 2aR 2 \og— Yu Dk (z) 

+a\z\ 2 + U»(z) 

= a\z\ 2 + aR 2 log r - P r 
Pi 
N' 
i? 

> F, (3 



F + 2aR 2 f 



where f{x) — x 2 1 — log x. Since / is nonnegative the effective potential satisfies 

U a (z) + V{z) > F. (3 



By Theorem 3.2 we conclude that the equilibrium measure for the background potential V is a. 
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□ 



The conclusion of Theorem |3.4| does not hold if some of the domains Dk overlap or intersect the 
exterior of B(0, R). In the first case we have to find a new decomposition of the perturbing measure and 
the corresponding domains; in the second the outer boundary no longer coincides with the boundary of 



B(0, R) as we saw in Prop. 3.3 It is hard to give a complete description of the support of the equilibrium 



measure in the general case. If v is a finite linear combination of point masses then the methods of 
D. Crowdy and J. Marshall in [3] used in the fluid dynamical context of rotating vortex patches are 
applicable to recover the corresponding supports. 

In the cases considered above, the support Sy of the equilibrium measure is contained the closed disk 
B(0,R). It seems plausible that the same is true for all perturbing measures v. The following theorem 



states that Sy c B(0, R) is valid if v is a positive rational linear combination of point masses. 
Theorem 3.5. Let V(z) be a potential of the form 



where v is a measure of the form 



V(z) = a\z\ 2 + U v (z) 



u = ^2 rk5a " 
fe=i 



(3.75) 



(3.76) 



where r%, r%, ■ ■ ■ , r m are positive rational numbers.. Then the support Sy of the corresponding equilibrium 



measure is contained in the closed disk B(0, R) where 



R 



l + i/(C) 



2a 



(3.77) 



Proof. For a continuous weight w(z) — exp(— V(z)) in the complex plane, z G C belongs to the support 
Sy if and only if for every neighborhood B of z there exists a weighted polynomial w n P n of degree 
degP„ < n, such that w n P n attains its maximum modulus only in B (see [15], Corollary IV. 1.4). Since 
our w(z) is continuous this characterization is applicable to this setting. 

Let z € Sy and suppose B is a neighborhood of z. Then there exists a polynomial of degree at most 
n for some n g N such that P n (z)w n (z) attains its maximum modulus only in B. Let q be the least 
common denominator of the rational numbers ri , r2, • • ■ , r n such that 



Pk 
<1 



(3.78) 



where p^eN for all k = 1, 2, . . . , m. Then 

{\p n { z )\w n {z)f = \p n ( z y\ 



Y[( z - a k) 



npk 



k=l 



-nqot\z\ 2 
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P n (z) q ~[[(z ~ a k ) n P* 



fe=l 
where 

m 

fe=l 

Since all the p/j's are assumed to be positive integers, 



exp ( ~n(q + L)^-\z\ 2 
q + L 



Q n{q+L) {z) := P n (z)i ~[[(z - a k y 



(3.79) 



(3.80) 



(3.81) 



k=l 



is a polynomial of degree at most n(q + L). If we consider the modified weight v(z) = exp (^— ^xl^| 2 ^ 
the corresponding weighted polynomial 

Q n(q+ L){z)v n{q+L \z) (3.82) 

attains its maximum modulus only in B. Therefore z belongs to the support of the equilibrium measure 
of the weight v(z) which is exactly B(0, R) where 



I? = 



I q + L _ / l + z/(< 
2qa ~ V 2a 



This proves that S v C B(0, R). 



(3.83) 



□ 



4 Orthogonal polynomials 

In random normal matrix models, the correlation functions are expressed in terms of planar orthogonal 
polynomials with respect to scaled weight functions of the form exp(— NQ{z)) associated to a potential 
Q(z) where N > is a scaling parameter (N has the same role as \). For our special potentials of the 



form V(z) — V a ^(z) defined in (2.14) above we have the weights 



e -NV(z) = exp (_ N r + r log 1 du ( w ) 

\ L J \z-w\ 

where N > is the scaling parameter. 
Proposition 4.1. We have 

e -NV(z) G L 1^ C) n L oo^ C) 

/or aZZ choices of the parameters N,a,v. Moreover, the absolute moments 

f \z\ k e' NV ^dm(z) 
Jc 

are all finite for k = 0, 1, . . . 



(4.1) 



(4.2) 



(4.3) 
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Proof. The exponent in the weight can be decomposed as 



N 



log- 



1 



\z — w 



dv{w) 



Na, 



\ 2 + N 



/log] — " — \ d v{ 
J \z-w\ 



w) 



in which the second term is lower scmicontinuous and satisfies 



N 



a, , 2 / , 1 



A + / log , 

2 ./ \z-w\ 



dv{w) 



N 



rv 1 

-|,| 2 + KC)log R 



O 



as \z\ — > oo. This means that the expression is bounded from below in C: 

1 



AT 



log 1 -dv{w) 

z — w\ 



> L 



for some constant L£l depending on the parameters N, a and on the measure v. So 

< e -NV{z) < e -^|,|^-L 5 

which implies that 



Moreover, 



for all k = 0,1,... 



/ |z| fe e- w ( z )dm(^) < e- L [ 
Jc Jc 



z | fe e N 2 a \ z \ 2 dm(z) < oo 



It follows from this and the positivity of the weight that the monic orthogonal polynomials 

P n , N {z) = z n + 0(z n - 1 ) (n = 0,l,-.-) 

satisfying 



Jc 



Pk,N{z)Pi,N{z)e 



-NV(z) 



dm(z) = h k ,N$ki, k, I = 0, 1, . . . 



(4.4) 
(4.5) 
(4.6) 

(4.7) 

(4.8) 
(4.9) 

□ 

(4.10) 
(4.11) 



exist and are unique where /i„,jv denotes the square of the L 2 -norm of P u ,n(z). 

5 Matrix d-problem for orthogonal polynomials 

In this section we show that the 2x2 matrix 9-problcm for orthogonal polynomials introduced by Its and 
Takhtajan [10] in the case of measures supported within a finite radius (cut-off exponentials of polynomial 
potentials) is also well-defined for the class of potentials considered above and determines the polynomials 
uniquely. In [10] the same family of potentials is considered as in [5]. 
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To be able to formulate the 9-problem, we need some estimates of Cauchy transforms of measures 



with unbounded supports. For a given potential V(z) — V aM {z) of the form (2.141 considered above let A 
be the measure that is absolutely continuous with respect to the Lebesgue measure in C having the form 

dX = e- NV ^dm (5.1) 

where and N > 0. Note that A(C) is finite because e~ NV W G L x {C,dm). 
The Cauchy transform of A is defined to be 

dX(w) 



[CX](z) := 



z — w 



(5.2) 



We need to control the asymptotic behaviour of the Cauchy transform [CA](z) of such measures at 
infinity allowing the possibility that [CA](z) is not holomorphic in any neighborhood of z = oo. 

that the density is bounded from above: e~ NV ^ < K for some 



4.1 



First of all, it follows from Prop. 
K G M. For a fixed positive radius R, we have 

dX(w) 



< 



< 



dX(w) 



1 



\z-w\>R 

dX(w) 



R 



Kdm(w) 



2tt r R 



\z-w\<R 



K 



-rdrdO 



A(C) + 1-kRK 



(5.3) 



for all z G C. Thus, there exists an upper bound H\ of [CA](z) depending only on iV, a, f and independent 
of z (One can get rid of R in the last expression e.g. by minimizing the bound in R): 

dX(w) 



Now 



[CX] (z) 



| z — w I 
A(C) 



< H x (z G C). 



(5.4) 



z — w 
wz 



dX(w) 



dX(w) 



1 



Hence the absolute value of the difference satisfies 

A(C) 



[CX] (; 



where A and A correspond to the measures 



< 



w\dX(w 
= I dX(w) + 



1 



z — w 
dX(w) 



dX(w). 



\w\ 2 dX{w) 



(5.5) 



< 



■Hx 



N 



S , 



(5.6) 



(5.7) 
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respectively. This means that 



[CX](z)= Xi ^+0 

z 



(5.8) 



If Pu,n{z) denotes the nth monic orthogonal polynomial with respect to the measure A, the modified 
measure 



d\ n {z) := \P ntN (z)\ 2 dX(z) 



corresponds to the perturbing measure 



(5.9) 



(5.10) 



k=l 



where {a^ 1 '^, a^"'^, ■ ■ ■ , «JT'' V '} are the zeroes of P„,jv(<z) and h n ^ = A„(C). An easy calculation gives 

1 



1 f»,w(^) - -Pn.Af(w) 1 Pn.Njw) 

z-w P n ,N(z) z-w Pu,n(z) Z-W 

-Q n - 1 (z,w) + 



(5.11) 



where Q n (z,w) is a symmetric polynomial in z and w of degree n — 1 with leading order z" _1 in the 
variable z. Therefore, by orthogonality, we get 



Ic 



Pu,n{w) 
z — w 



dX(w) 



— I 

°n,N(z) Jc 



\Pu,nH\ : 

C z - w 
dX n (w) 



dX(w) 



Pu,n{z) J c z-w 



1 z n 
z"P n , N ( z ) 

7 n+l 



h n . 



N 



+ o 



y n+2 



O 



(5.12) 



Following the approach of Its and Takhtajan in [10], we consider the following 2x2 matrix- valued 
function in the complex plane: 



Y k<N (z) := 



Pk,N(z) 



Pk ^e- NV ^dm(w) 
ft Jc w — z 



lk-l,N tlk-l,N Jc W — Z 



m(w) 



(5.13) 



The 9-derivative is 



^Y KN {z) 



-P k , N (z)e 




hk-i,N 



NV(z) 

P k . ltN (z)e- NV ^ 



Yk,N{z) 



-e~ NV W 




(5.14) 
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Using the asymptotic behaviour of the Cauchy transforms as \z\ — > oo proven above, we have that 
Pk, N (z) 



1 f P^H e -NVM dm{w) 



l k-l,N 



-Pk-l,N{z) 
Pk,N{z) 



l k-l,N 



71" Pk-l,N(z) 



z 



hk-i.N 
1 



hk-l,N Jc w 



Pk-l,NH e -NV {w)dm{w) 

w — z 

w — z 

Pk - lMw) e- NV ^dm(w) 











1 



+ o 



(5.15) 



So Yfc,Ar(z) is a solution of the following 2x2 matrix d- problem: 



— M(z) = M{z) 
oz 



M{z) 



-NV(z) 





" 1 


" 


( 





1 




z k 
z~ k 



(z E C) 



(N^oo). 



(5.16) 



The important point made in [10] is that the 9-problem in that setting has a unique solution and 
therefore it characterizes the matrix Yk^(z) and the corresponding orthogonal polynomials. Although we 
cannot assume that the relevant Cauchy transform entries are holomorphic around z = oo, we nevertheless 
can prove that the solution is unique in this case as well. 



Proposition 5.1. The matrix Y^^(z) is the unique solution of the d-problem (5.161. 



Proof. We have seen that Yk,N solves the <9-problem (5.161. Conversely, assume that the matrix M(z) 



has continuous entries with continuous partial derivatives and M{z) solves (5.161 with the prescribed 
asymptotic conditions. Then Mu(z) and Mi\{z) are entire functions with asymptotic forms for large z 



M u (z) = z k + 0(z k - 1 ), 
M 2l {z) = 0{z k ~ 1 ) \z\ 



oo. 



(5.17) 
(5.18) 



Hence Mn{z) is a monic polynomial of degree k and M 2 i(z) is a polynomial of degree at most k — 1 

en in tc 

M 12 (z) 



The 9-equation in (5.16 1 can be written in terms of the entries of M(z) as 

d_ 



d_ 

dz 



M 22 (z) 



(5.19) 
(5.20) 
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Taking into account the fact that Mi 2 (z) — > and M 22 (z) —* as |z| — > oo this implies 

Ml2 (z) = 1 / ^lM e -w( w ) dro(lu)j (5 21) 

M 22 (z) = 1 / ^iM e -^W dm(w) (5 . 22) 



(see [2]). Using the expansion 



1 1 z k — w k 1 w k 
+ 



z — w 

fc-1 



we get 



Ml2 (,) = I / ^ e - w Mdm(») 
n J c w- z 

k — 1 

Y 4rr- f w l M^fw)e- NV ^dm(w) 



i=0 



11 r u, MnW e _ NV{w) 

Z K IT J c W — Z 

The prescribed asymptotic behaviour 

M 12 (z) - O (^tt) as |z| - oo (5.25) 



implies the following equations: 

/ w / M 11 (w)e- Ary(tu) rfm(w) = Z = 0,1,..., fc-1. 

Hence Mn(z) = Pk,N{z) because Mn(z) is a monic polynomial of degree k. Similarly for M 2 i(z) we 
have 

M 22 (z) = 1 / M^R e -NV iw)dm{w) 

IT J c W-Z 
k — l 

= V^tt- I w l M 21 (w)e~ NV ^dm(w) 



11 f ^W e^) (5 . 26) 



Z k IT J r W — Z 



and 



implies 



M 22 (,) = l + o(-l T ) as|*|-oc (5.27) 
/ io'M 2 i(«))e- w W(im(w) = J = 0,1, . . . , A - 2, (5.28) 
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and 

' w^M^Jwje-^^dmiw) = 1. (5.29) 



c 



Now, if M 21 (z) = az k - x + O (z fc ~ 2 ), where a £ C, then 

\M 21 (w)\ 2 e- NV ^dm(w) = a [ w^ 1 M 21 (w)e- NV(w Um(w) = a. (5.30) 



Clearly a ^ (because otherwise M2i(z) and hence also M22(z) = would be zero which is impossible). 



So M 2 i(z) is a polynomial of degree k — 1 and from (5.281 we have that 



Af 21 («) =aP fc _i, JV (z). (5.31) 



By the asymptotic relation (5.121 



M 22 (z) = - [ ^lM e -NVM dm{w) 



7T <n z — w 



ahk-i.N 1 , m ( 1 



^ ( Zk+i ) ' (5-32) 



which forces the constant a to be equal to — j— 1 — and hence 



which completes the proof. 



M 21 (z) = P fe _ 1JV (5.33) 

ftfc-l,JV 



□ 



6 Zeros of Orthogonal polynomials and quadrature Domains 

In this final section we briefly discuss some known and conjectured relations between the asymptotics 
of orthogonal polynomials in the plane, equilibrium measures of the type studied in Sections 1-3 
and generalized quadrature domains. These concern relations between the asymptotics of the zeros of 
orthogonal polynomials and the associated equilibrium measures that have previously been studied by 
Elbau [5, 6] for a certain class measures with bounded support. Their validity for the class of measures 
considered here is supported by numerical calculations. 

To relate orthogonal polynomials with the support of the equilibrium measure we have, for finite 
values of n, three quantities which, in the case of random Hermitian matrices are known to approach the 



same equilibrium measure on the real line in the scaling limit (cf. ( 1.3 1 ) 



oo, 7V->oo, —^ 7 =i, /j := _L e M+. (6.1) 
n't N v ' 
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All the limiting relations below are understood in this sealed sense. We introduce the notation 

Q(z) ■■= \v{z) (6.2) 

for the rescaled potential corresponding to the scaling parameter 7. Then for a large class of real potentials 

and z e R, all three of the following measures converge weakly to the equilibrium measure dfiQ(z) of 
Q(z): 

1) The normalized counting measure of the zeros 

Z n , N :={z[ n ' N \zi n ' N \...,z^} (6.3) 
of the orthogonal polynomials P n ,N(z) with respect to the weight exp(— NV(z)) 

v n ,N ■=- E 5 Z , (6.4) 

v n ,N (n-»oo). (6.5) 

2) The expected density of eigenvalues (or one-point function) of random Hermitian matrices 



1 " * 

P»A*) = " E \PkM\ 2 e~ NV(z \ (6.6) 

fc=0 

drawn from the probability density 

J—eM-NTv(V(H)))dH, (6.7) 

Z n M := J cxp(-NTr(V(H)))dH (6.8) 

p n ,N(z)dz dfi Q (z) (n -> 00). (6.9) 
3) The normalized counting measure of equilibrium point configurations 

J*:={z?,...,z*} (6.10) 

of the two-dimensional Coulomb energy 

N N 

£ n (z u ... ZN ) = - io gl —— ,+eq(*o ( 6 - n ) 

i,i=i > Zt z i\ i=l 
(which in the plane become the so-called weighted Fekete points) 



Vr 

n 



I E *. ( 6 - 12 ) 
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Vn — > MQ (n— ►<»). (6.13) 
For random normal matrices the eigenvalues are not confined to the real axis. In this case it is 



known [1,6,9] that in the scaling limit (6.1 1, the analogs of p n ^{z)dz and rj n also approach the equilibrium 
measure \iq. 

It is also known [6] that for cut-off measures of the form 

e- NV ^I v {z), y(z) = - a |z| 2 + P harm (z), (6.14) 

where I-p is the indicator function of a compact domain T> containing the origin whose boundary curve 
dT> is twice continuously differentiable and Pharm is a harmonic polynomial, that 

1 1 f 1 
Jim -log — — - = / log, -g^q(C), zeC\supp(/j Q ). (6.15) 

That is, the limit of the zeros acts effectively as an equivalent source of the external Coulomb potential. 

For such potentials, the boundary <9(supp(/iQ)) of the support of the equilibrium measure is determined 
through the Riemann mapping theorem as the image of the unit circle under a rational conformal map, 
whose inverse therefore has a finite number of branch points. The Schwarz function S(z), defined along 
the boundary, determines the curve via the equation 

z = S{z). (6.16) 

It has a unique analytic continuation to the interior on the complement of any tree C trec whose vertices 
include the branch points, with edges formed form curve segments. It is shown in [6] that, assuming 
there is a condensation limit for the orthogonal polynomial zeros supported on a tree-like graph Cq 
whose edges are curve segments, C troG may be chosen so that C C troc . Moreover, denoting by SS(z) 
the jump discontinuity of S(z) away from the nodes, C trcc may be chosen as an integral curve of the 
direction field annihilated by the real part Ke[(SS(z)dz] of the differential (5S(z))dz; i.e. such that the 
tangents X to the curve segments forming the edges satisfy 

Re[{6S(z)dz](X) = 0. (6.17) 

We refer to such integral curves as critical trajectories. 

Based on computational evidence, and general results known for other cases [13], there is good reason 
to believe that the same result holds for the class of superharmonic perturbed Gaussian measures studied 
in Sections 1-3, without the need for introducing the cutoff factor /-p. This statement, for some suit- 
able restrictions on the permissible superharmonic perturbations, forms the first part of the conjectured 
relation between the zeros of the orthogonal polynomials considered in Sections 4-5 and the equilibrium 
measure hq. 
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The second part gives a more detailed relation; namely, the effective density Kq(z) along Cq of the 
condensed orthogonal polynomial zeros is given, within a suitable scaling constant by 



dK Q (z) ~ ±(SS(z))dz = ^lm[(5S(z))dz}. 



(6.18) 



Explicitly, this means that the external potential due to a uniform, normalized charge supported in 
supp(a*q) is 

log ^—.dngiC) = [ log - J—d^fl. (6.19) 

Cq N-CI J I 2 "CI 

To support the validity of these conjectures, we take the case of the potential 

1 



V{z) =a\zY + p\o£ 



(6.20) 



\z — a\ 

in the simply connected case considered in Section 3 above and compare the locus of the zeros of the 
corresponding orthogonal polynomial P n n{z) with the two different integral curves of the direction field 



defined by (6.17) joining the branch points (Fig. 6.1). (The two other critical trajectories emanating 



from the branch points are omitted from the graph) 




Figure 6.1: Zeroes of P n ,N(z) for n = 50 and the critical trajectories 



We also compare, in Fig. 6.2, the value of the logarithmic potential — - log |P n ,Jv(2)| created by the 
normalized counting measure v n ,N of the zeroes of P n ,N(z) with n — 30, N = In — 60 and the external 



potential as given by (6.19) in the external region zeC \ supp(//y). 
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-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 -l'.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

x x 

Figure 6.2: The figure on the left is a contour plot of the potential — - logP ni Ar(2:)| arising from 
equal charges at the roots of the orthogonal polynomial P U} n(z); the one on the right is the poten- 
tial /log d/iQ^) on the exterior region due to a normalized uniform charge on supp(/zq). 
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